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\ We calculate the effect of an inhomogeneous magnetic field, which is invariably present in an 

■ experimental environment, on the exchange energy of a double quantum dot artificial molecule, 

' projected to be used as a 2-qubit quantum gate in the proposed quantum dot quantum computer. 

We use two different theoretical methods to calculate the Hilbert space structure in the presence 
of the inhomogeneous field: the Heitler-London method which is carried out analytically and the 
molecular orbital method which is done computationally. Within these approximations we show that 
the exchange energy J changes slowly when the coupled dots are subject to a magnetic field with a 
wide range of inhomogeneity, suggesting swap operations can be performed in such an environment 
as long as quantum error correction is applied to account for the Zeeman term. We also point out 
the quantum interference nature of this slow variation in exchange. 

'd 

r-j ; PACS numbers: 03.67.Lx, 73. 21. La 



(N 



C/3 



S 

I 



> 

o 



I 



- 

X 



I. INTRODUCTION 



It has long been proposed that quantum mechanical computers might provide much more computing power than 
machines based on classical physics [Q. However, it was the discovery of fast quantum algorithms for tasks such as 
■ prime factorization and searching disordered databases [|| together with effective error correction schemes |^ that 
Ch ' motivated a frenetic search for a suitable physical system on which reliable quantum hardware could be built. 
O , Among the many proposed quantum hardware schemes the ion trap quantum computer, based on electronic states 
I ! '• of laser cooled trapped ions, have not only demonstrated full control of elementary one and two-qubit (quantum bit, 
the basic unit of quantum information) logic gates, but also recently produced 4-qubit entanglement Q. Although the 
ion trap represents the current state of the art in quantum computing architecture, its scalability prospect is unclear, 
since moving a large number of ions between memory and computing sectors of the ion traps without destroying 
quantum coherence is a difficult problem. Another proposed scheme, the liquid state Nuclear Magnetic Resonance 
' (NMR) quantum computer, based on nuclear spin states of molecules in solutions, recently demonstrated a 5-qubit 
order finding algorithm However, the NMR quantum computer also has a scalability problem, since the signal to 
noise ratio decreases exponentially as the size of the molecules increases. Furthermore, there has been a theoretical 
T-H proof that no entanglement is present in the current ensemble-averaged NMR experiments Q , although its implication 
for NMR quantum computation has not yet been fully explored. 

With the scalability problem in mind it has been suggested that advanced semiconductor technology might provide 
the basis for building a large (i.e. many qubits) quantum computer (QC), just as it has been proved to be effective 
in building today's classical digital computers. Two of the main proposals for such semiconductor QCs are based on 
qubits using the spin of a phosphorous donor nucleus in bulk silicon |^] and the spin of a single electron trapped in a 
^ • GaAs quantum dot ||^. Even though our work here refers specifically to the latter, our results are valid for both solid 
Q , state QC proposals since the two are similar in the sense that they use the same universal quantum gate to perform 
O ■ quantum computation: The exchange or swap gate which acts on two qubits coupled by an effective Heisenberg spin 
interaction and one-qubit gates that perform rotations on two orthogonal spin axes by applying external magnetic 
fields. Experimentally demonstrating the operation of any of these quantum gates is a formidable task and has not 
been carried out yet, and all semiconductor based QCs are at best promising proposals at this stage. Nevertheless, the 
hope is that after the first step of successfully demonstrating a single qubit in a semiconductor-based QC, scaling up 
from one qubit to many would be feasible based on the established infrastructure of semiconductor microelectronics 
industry. Moreover, precisely controlling the quantum dynamics of a few spins in a semiconductor nanostructure 
represents a major experimental challenge, and therefore any success will lead to new insight into mesoscopic spin 
properties of nanostructures, which by itself is a subject of intrinsic fundamental interest. 

Quantum computation can be considered as a complicated unitary evolution of a N-qubit system, governed by 
a N-qubit Hamiltonian, where each qubit is a quantum mechanical two-level system. It has been proved that this 
usually extremely complicated process can be factored into a number (that is a polynomial in N) of simple one and 
two-qubit processes with sufficiently high accuracy, in analogy with a classical network of logic gates H . To improve 
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efSciency multiple simple operations might be lumped together, as was proposed in various parallel pulse schemes 
1^^ . In these operations for a quantum dot quantum computer (QDQC) an inhomogeneous magnetic field inevitably 
appears, as it is a powerful technique to rotate spins while controlling spin exchange. The inhomogeneity in the field 
here arises necessarily from the fact that the parallel scheme employs local magnetic field pulses on individual dots 
to achieve different spin rotations, and a magnetic field localized on one dot is essentially an inhomogeneous field 
for a different dot. Furthermore, inhomogeneous magnetic fields can be produced by unwanted magnetic impurities 
or external currents in the sample, which will be invariably present in any real QDQC architecture. Since in the 
proposed parallel schemes the exchange coupling J (which controls the swap gate) is treated as independent of the 
magnetic field inhomogeneity pO| , it is important to explore whether this dependence is truly weak enough, or if it is 
necessary to adjust the inter-dot barrier continuously in order to maintain the desired J in the presence of a varying 
inhomogeneous field (which will be a formidable, if not impossible, task). 

We showed in a previous paper that an inhomogeneous magnetic field leads to errors in the swap operation of 
two qubits because the Zeeman term in the effective Hamiltonian is not proportional to the total spin. For example, 
if we try to swap the state |ti) we would obtain ||t) plus an entangled component proportional to the magnetic field 
inhomogeneity. In other words, once a product state is entered, the output will always be entangled, although for 
practical purposes this error can be kept within the currently estimated error correction bounds (~ 10^^) for usual 
field inhomogeneities. This result is based on the assumption that one knows J exactly in order to achieve the best 
control, otherwise there can be additional errors in swap because of the uncertainty in J. 

Motivated by the above concerns and problems arising from inhomogeneous magnetic fields, we calculate in this 
paper the exchange coupling J of a two-electron double dot as a function of the inhomogeneity in the external 
magnetic field. We employ a single envelope function approach and calculate the orbital energies of the two electrons, 
thus relating the general two-electron problem to a special two-spin Hamiltonian, and analyze the effect of field 
inhomogeneities on the orbital degrees of freedom of the double dot structure. The current work is thus complementary 
to our recent publication dealing with the field inhomogeneity induced errors in the swap operation. 

Quite recently it has been demonstrated that QDQC can be done without any applied magnetic field as long as a 
price is paid in increasing the number of exchange operations: Two {S = ^, Sz = +5) states of three spins are used as 
a qubit, and the nearest-neighbor exchange interaction can provide all the required one-qubit and two-qubit operations 
of a quantum computer fl^ . In this scheme applications of external magnetic fields for single qubit operations is not 
required, as exchange gates provide both 1- and 2-qubit operations. This proposal removes the problem of having 
strong applied inhomogeneous fields. However, stray fields due to impurities and external currents can still be present, 
and must be taken into account using appropriate error corrections. In this paper we focus on clarifying the issues that 
arise in the simpler sequential scheme that may contain an applied inhomogeneous magnetic field as well as pointing 
out potential errors in all exchange QC schemes. Our results are also directly relevant to the proposed parallel pulse 
schemes as discussed above. 



We consider two electrons in a horizontally coupled GaAs double quantum dot structure (Fig- !)• Since 

these electrons are at the bottom of the conduction band the single envelope function approach provides a sufficiently 
complete physical picture of their energy states. Suppose there is an inhomogeneous magnetic field applied so that the 
electron in the left well is subject to a field which is on average 2Bdif lower than the one in the right. For simplicity 
we assume a linearly changing field profile 



II. MODEL 




where B — ^ {B^ + Bl), Bdi f — ^ [Br — Bl ) , and 2am is the separation between the dots . The complete Hamiltonian 
is given by 
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Hz = g*iiB ^ B(xi) • Si. 



(6) 



Here is the position of the ith electron in the xy plane, m = O.OGTtoe is the electron effective mass in GaAs, 
e = 13.1 is the dielectric constant, g* = —0.44 is the GaAs gyromagnetic ratio, and A is the vector potential for the 
inhomogeneous field: 
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Note that since the quantum dots are fabricated using a high mobility 2D electron gas (taken to be along the xy 
plane) , the dot dynamics is assumed to be constrained to the z = plane because the z- width of the 2D electron gas 
is typically much smaller than the QD confinement in the x-y plane. The z-dependent components in A and B do 
not have any effect other than to guarantee the V x B = condition. Hence the effect of Bx = Bdifz/um will be 
neglected in our calculation (the 2D electron gas has a typical thickness of the order of a few nm and neglecting B^ 
is an excellent approximation for realistic QD systems). The confining potential V in the x-y plane is modeled in our 
calculations by a linear combination of gaussians characterized by the parameters Vq, a, Vf,, l^, ly, hx and 4y which 
can be adjusted to control the size of the dots and switch on and off the exchange coupling: 
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In our calculation we choose Vq = BOmeV as the potential well depth and 14 = SOmey as the central barrier height. 
In addition, a — ISnm, Ix — ly — 30nm, Ibx = ISnm, and 4y = SOnw, which leads to an effective barrier height 
(difference in energy between the top of the central barrier and the bottom of one of the wells) of 9.65meV, and 
a slightly shifted potential minimum a„i — 16.417ito. We note that our choice of Eq. (||) as the QD confinement 
potential model is actually quite a general approximation, and allows considerable flexibility in our QDQC Hilbert 
space structure calculations. 

An important feature of the Zeeman Hamiltonian is that it couples coordinate Xi with spin Si in an inhomoge- 
neous field. This interaction introduces a coupling between the singlet and the 5"^ = triplet state, and can be well 
described (in the Heitler-London approximation) by an effective spin Hamiltonian |ll|] 



Hz = 7rSiz + 7lS2z, IR = g*fJ,BBR, 7l = g*l-J.BBL- 



(9) 



But note that Eq. (|^) is not invariant under particle permutation P12 , which is a natural consequence of the fact that 
it has a non- vanishing matrix element between singlet and triplet states (which are respectively odd and even under 
a P12 operation acting only in the two-spin space) . This lifting of the permutation symmetry is a direct consequence 
of the inhomogeneous external magnetic field in this problem. 

Our main objective is to study the dependence of the exchange coupling J on the inhomogeneity of the magnetic 
field. The parameter J describes the splitting between the lowest energy singlet and triplet states of the Hamiltonian 
(^, which at low energy reduces to the effective two-spin Heisenberg-type Hamiltonian that enables the operation of 
the quantum gates. 



Heff = J Si • S2 + ^rSiz + JlS2 



(10) 



We use four different approximations to calculate J from the orbital Hamiltonian given in Eq. (j^): Heitler-London 
(HL), Molecular Orbital with S states (MO-S), with up to P states (MO-SP), and with up to D states (MO-SPD). It 
should be noted that the Zeeman splitting term given in Eq. (^ does not affect the calculation of J in any way and 
it is only added to Eq. (10) after J has been determined in order to obtain the effective 2-spin exchange Hamiltonian 
controlling the quantum gates. 



III. METHODS AND RESULTS 



To diagonalize the Hamiltonian given in Eq. (^, we use as our one particle basis the so called Fock-Darwin energy 
states centered at the bottom of the left and right wells {x = ±a„i, y = 0). The Fock-Darwin QD states, briefly 
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described in Appendix A of this paper for the sake of completeness and notational clarity, are simply the quantized 
states of the 2D orbital motion of the electrons in the presence of an external parabolic confinement potential and 
magnetic field. Fock-Darwin states have a principal quantum number n which describes a shell structure: n=0 
is referred to as S states, n=l as P, n—2 as D, and so on. Each shell is composed of n+1 levels labeled by the 
magnetic quantum number m=-n, -n+2, n-2, n. Their energy is given by Enlrt = TiVt^/^(n + 1) + ^hiOc'^^m where 

^R/L _ Y^|^2 _|_ {ijj^/^ /2Y^ ijJc^^ = eBjiii^/mc is the cyclotron frequency and luq = y^Vo/ml'^ approximates the 
harmonic confinement frequency near the bottom of the dot potential. Each Fock-Darwin state has a characteristic 
width Iqj^/l — y'hJrnSW/^ , which assumes different values depending on whether the state is in the right well subject 
to a field Bn or in the left well where the field is B^, which are different in the inhomogeneous field case we are 
studying. Notice that the vector potential (0) is centered in between the dots which is located ±0^ to the side of 
each well. This means that we have to use basis functions which are gauge transformed Fock-Darwin states with an 

additional magnetic phase exp ^±z— ^^^^^^ (sign -I- and — refer to the left and right wells respectively). 

In what follows we use this Fock-Darwin single electron basis to form two-electron states and solve our two-particle 
problem in two approximations: Heitler-London which can be solved analytically and Molecular Orbital which is 
solved numerically. 



A. Heitler-London 



Heitler-London (HL) is the simplest method to calculate the exchange coupling J in two coupled dots. It consists 
of the singlet and triplet wave functions formed from the lowest Fock-Darwin states (n=m=0): 
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where the spin functions are Xs = ^ (|Ti) - lit)) and Xt is one of the three triplet states ^ (|ti) + lit)) , ITT) , lii)- 
The exchange energy is then calculated as the difference between the expectation values of the Hamiltonian (|^) 
calculated with the singlet and triplet trial wave functions: 



J ={T\Hort,\T)-{S\H,rb\S) 



(13) 



Using this formula in an inhomogeneous field, we can obtain an exact expression for J, which is quite complicated 
when both B and Bdif are nonzero, and is given in its entirety in Appendix B. Here we show the expression for J in 
the special situation when B — Q, which catches the main qualitative features of the inhomogeneous case: 
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where S = exp (— am/^o) is the overlap, Iq 



Iqr — loL the Fock-Darwin characteristic width (the same for the left 
and right wells in the B = case as \Bl\ = |i?fl|), and Iq the zeroth order modified Bessel function. Note that in 
general am 7^ ci since the location of the minimum in the wells is usually different from the parameter a = 15nm. 
The main difference between Eq. (|lj) and the homogeneous field Bdif = 0, B ^ case (see Appendix B) is that the 
latter contains a number of phase interference terms proportional to exp [— (ei3/2?ic)^a^ZQ] , responsible for the faster 
decay of J as a function of B as will be discussed below. In Fig. 2 we depict, primarily for the purpose of comparison 
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with the inhomogeneous field results shown in Figs. 3 and 4, the calculated exchange coupling J as a function of the 
external magnetic field B in the homogeneous field case with B^if = in all our approximations. 

Two curves calculated for different inhomogeneities using the HL approximation are shown in Fig. 3. The HL 
method is qualitatively reasonable, however, quantitatively it differs significantly from the more sophisticated molec- 
ular orbital calculations described below in section III B because J depends sensitively on the wave function overlap, 
which can be significantly affected by including higher excited states (as is done in MO calculations, but neglected in 
the HL theory). Even though the maximum field gradient obtainable in a laboratory is around 1 T in a 30nm range, 
we show values of Bdi / that go all the way up to 5 T for completeness. 



B. Molecular orbital calculation 



The basic idea of the MO calculation is to mix in higher excited states in the set of basis states so as to improve 
the energetics. To go beyond the HL approximation, we can assemble more two-particle wave functions from the 
Fock-Darwin states, which we call molecular orbitals 

vl/,(l,2) = (*,(i)(l)vI/,(2)(2) ± M/,(2)(l)vI/,(i)(2)) Xs/t . (15) 

Electrons are fermions so the singlet has to be symmetric in space while the triplet is antisymmetric. Here we use 
this molecular orbital (MO) basis to solve a generalized eigenvalue problem H • Vl/ = S • \I' where Hij = \ H 

= (5*1, ^"2, . . .) is the vector state and the overlap matrix Sij =< 'i'il'^j > is needed for orthonormalization since 
the MO basis is in general not orthonormal. We have done such a calculation in a constant homogeneous field p4[ , 
using both S and P Fock-Darwin states. Here we perform calculations with the progressively larger basis sets of S, 
SP, and SPD states in the inhomogeneous field problem. 

In the simplest MO calculation we use only the Fock-Darwin S states to form the two-electron orbitals. Such a 
calculation improves upon the HL method by incorporating the doubly occupied states ^_R(l)\l/ij(2) and 5'l(1)^'l(2) 
and we call this the MO-S calculation. We then have a 3 x 3 singlet matrix to diagonalize, and the triplet energy is 
the same as the HL case (a 1 x 1 matrix). Note that we cannot use parity symmetry here to break up the matrices in 
a block diagonal form as was done before since the field is inhomogeneous in the x direction. The resulting MO 
calculation is carried out numerically. The results of such calculations are shown in Figures 2 (homogeneous case), 
and 3 (inhomogeneous case) where there is a noticeable quantitative difference between the HL and MO-S exchange 
energies. 

To improve the accuracy of our molecular orbital calculations, we also perform an MO-SP calculation, adding to 
our basis the n = 1, m = ±1 states on each side so that we have 6 one particle levels in total. The singlet MO-SP 
/ 6 \ 

matrix is 21x21 (( ^ 1+6 = 21) and the triplet matrix is 15 x 15. We are able to efficiently calculate the Coulomb 

matrix elements Cy — (^'^j C j'l'j) by reducing the 4-dimensional two-particle integrals to a ID integral using the 
identity 

■.expi-srl^) (16) 



|xi-X2| ri2 \^ Jo ^/s 

and integrating over Xi and X2 in the gaussian basis ||l^. We then further increase the size of the MO basis by 
including Fock-Darwin D states (n=2, m=-2,0,2). We now have 12 single-particle levels, leading to a 78 x 78 and a 
66 X 66 two-electron Hamiltonian matrix in the singlet and triplet basis, respectively. 

Before we analyze our results it is important to determine the range of validity of our approximations. A strong 
inhomogeneity across the region where our basis functions vary appreciably leads to the necessity of including more 
and more orbitals to calculate the energy levels accurately. The extra terms that a non-vanishing inhomogeneity 
{Bdif 7^ 0) adds to the one particle Hamiltonian can be written explicitly as a one particle operator AH. We write 
the vector potential as A = Ar + Sh where Sj^ represents the deviation from An = —^yBBX+ ^xBjiy which describes 
the field in the right well. The one particle Hamiltonian becomes 

h^ = :^(P. + -Ajif + V + —SR.ip, + -An) + :^i-fSl . (17) 
Zm c mc c Zm c 



AH 



AH is the part of the orbital Hamiltonian (m) that is most sensitive to inhomogeneities in the magnetic field, since it 



is proportional to B^^r (the Coulomb matrix elements do not depend strongly on Bdif). In an inhomogeneous field. 
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Aff strongly couples the S orbitals with the P and D orbitals, as can be seen from Table 1, where this orbital mixing 
is quantitatively represented. Each individual term in A_ff contributes approximately the same to the coupling, so we 
cannot neglect any of the terms in the calculation. 

Note that as Bdif increases, the corrections to the MO energy levels grows very fast, reaching 10% of the exchange 
coupling J ( J = 0.2meV a,t B = OT and J = O.lmeV at B — 4T) for B^if = IT in the second and third columns of 
Table 1. This suggests that the HL, MO-S and MO-SP approximations are only reliable up to B^if = IT for both 
the B = OT and B = AT cases. Moreover, a good estimate for the reliability of our best approximation, MO-SPD 
can be set as Bdif = 4T for the B = OT case and Bdif = 2T for the B — AT case. In Fig. 3 we show the exchange 
coupling J calculated only up to Bdif — 5T, which is an upper limit on the reliability of all our approximations. For 
completeness. Fig. 4 displays the same J calculated up to Bdif — 25T (In this context, we should emphasize that in 
realistic experimental situations Bdif is likely to be much less than 1 T). 

On the other hand, our approximations are increasingly better at low inhomogeneity, as can be seen in Table I. 
Fig. 3 shows that the exchange coupling J decreases only by 0.1% in the B — OT case when Bdif = 0.5T, which 
is the maximum field gradient obtainable in an experiment {dBz/dx = 300Gauss/nm). Also, in the B = AT case 
it increases by 0.6% when the same field gradient is present. Using these values we can estimate that the error 
correction upper bound (10~^) corresponds to an inhomogeneity of 30Gauss/nm. In section IV below we will argue 
that quantum interference in the inhomogeneous case is responsible for such a small change in J, as contrasted with 
the homogeneous situation in Fig. 2. Moreover, Fig. 3b shows that the exchange coupling J slightly increases as we 
switch on the inhomogeneity in a finite field {B = AT). All our approximations yield an increasing J, which can be 
understood in a simple one-particle picture. As we increase Bdif the magnetic field on the right (-B/j) also increases, 
squeezing the Fock-Darwin state located on the right quantum dot. However, the left magnetic field (B^) decreases 
to less than 4 T, which leads to a broader electronic wave function on the left quantum dot. When Bdif = 4T, the 
electronic state on the left dot reaches its maximum width, and starts to be squeezed for greater inhomogeneities. 
This suggests a maximum of the exchange coupling J located around Bdif = 4T. This maximum is not evident in 
our approximations simply because they lose their accuracy for Bdif below 4 T as discussed above, but an increasing 
trend in J is apparent from Fig. 3b. Moreover, by looking at the B = OT case (Fig. 3a) we can infer that J is probably 
not increasing for Bdif < 2T, since, as the MO approximations become more accurate, J first increase (MO-S and 
MO-P) and then decreases (MO-SPD). No maximum in J is expected here since as Bdif increases both the left and 
right electrons are squeezed. 



IV. DISCUSSION 



Comparing the results presented in Figs. 2 and 3 it becomes evident that J decreases much faster (Fig. 2) in the 
homogeneous magnetic field, even becoming negative at around 5T, in both HL and MO approximation schemes. 
This behavior can be roughly understood as arising from the approximate proportionality of J to the wave function 
overlap S = J ^'^^E'^dr {'^r/l are defined in Eq. (p^) multiphed by a polynomial in B and Bdif which explains 
features such as the singlet-triplet crossing in Fig. 2 (the overlap S is plotted in dimensionless units together with 

our results for J). It turns out that S is very sensitive to the magnetic phase exp (^zLi—^^^^^ which appears 

in the right and left basis functions discussed above. In a homogeneous field {Bn — Bl = B) those phases add 
up, leading to destructive interference which gives an extra exponential decay proportional to B^ in the overlap: 

SHom = exp (— (flrn/'o)^) ^xp ^— (s^o'^mf^)^^ ■ On the other hand in the inhomogeneous case with 5 = 0, Br = —Bl 

the phases cancel out in the S integral leaving just a single exponential decay. Sink = exp (— (am/Zo)^) which depends 
mainly on the interdot distance. This absence of quantum interference in the inhomogeneous field case explains why 
J hardly changes when there is a magnetic field difference over the two quantum dots (Fig. 3), which is in fact good 
news for the operation of the swap gate in an environment that contains some magnetic impurities for example (or a 
field inhomogeneity arising from some other unintended source). 

Note that in the homogeneous case the exchange coupling constant in Fig. 2 decreases extremely strongly at a 
few Tesla which leads to a singlet-triplet crossing. This crossing arises from the competition between the long range 
Coulomb interaction (which favors the triplet state) and the kinetic energy (favoring the singlet state). There is no 
such singlet-triplet crossing in an inhomogeneous field (except when B is such that J < 0, where the application of 
Bdif leads to a positive J), because the kinetic energy term (which favors the singlet state) now dominates as quantum 
interference terms are suppressed. In the homogeneous field case on the other hand the potential energy term always 
wins at high enough B fields leading to the triplet state becoming lower in energy than the singlet state. 
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V. CONCLUSION 



In conclusion, we study the singlet-triplet splitting (i.e. the exchange coupling) J of a horizontally coupled double 
quantum dot system with two electrons. In particular, we calculate J when the magnetic field is inhomogeneous with 
different field strengths {AB = 2Bdif ) in each dot, and show that realistic values of field inhomogeneity in double 
quantum dot structures should have quantitatively small (but qualitatively interesting) effect on J. We perform 
an analytical Heitler-London calculation, and a numerical diagonalization of the orbital Hamiltonian in a set of 
increasingly larger two-electron MO basis to conclude that the cancellation of the magnetic phase difference between 
the two electrons on the dots is responsible for the relatively weak dependence of J on Edif for realistic values of 
Bdif. The swap gate can thus be operated in an inhomogeneous environment as long as the exchange coupling J is 
properly chosen and quantum error correction is performed to account for the inhomogeneity-induced errors arising 
from the Zeeman term as discussed earlier We conclude therefore that electron spin based QDQC operations 

should in principle be feasible in an inhomogeneous magnetic field environment (with the inhomogeneity arising either 
from external magnetic fields needed for single qubit manipulation, or for running the parallel pulse scheme, or from 
unintended magnetic impurities or currents invariably present in the system) as long as the field inhomogeneity is not 
too large. 

This work is supported by ARDA and the US-ONR. RdS acknowledges useful discussions with Dr. J. Fabian. 



APPENDIX A: FOCK-DARWIN STATES 

We briefiy describe the Fock-Darwin states which are used as the single particle basis in our calculations. 
Consider an electron of mass m* restricted to move in the x-y plane and confined in a parabolic well of frequency wq. 
A magnetic field B is applied in the z direction, which can be described by the vector potential A = — i?/i3x-|- ^xBy. 
The Hamiltonian then becomes 

where r — \/ + is the distance between the electron and the origin, e the electronic charge and c the speed of light. 
Since this Hamiltonian commutes with the z component of the angular momentum operator Lz, its eigenfunctions 
can be chosen with well defined angular momentum mh. Apart from the magnetic quantum number m, a radial 
quantum number n completely describes the eigenenergies of the system, given by Enm — {n + l)hil + ^mhojc, with 
71 = 0, 1, 2, . . . and m = — n, —n + 2, . . . , n — 2, ri. Here ujc — is the electron's cyclotron frequency whereas 
i7 = -y/wg -f {ijJc/2)'^ is the harmonic oscillator frequency renormalized by the magnetic field. For small magnetic 
fields {u)c "^o) the quantum number n describes a shell structure, with n = 0, 1, 2, . . . being the S,P,D, . . . shells 
respectively. In table 2 below we write the angular momentum and energy eigenfunctions for the S, P and D shells 
that are used in calculations throughout this paper. The negative m wave functions can be easily obtained from the 
positive ones using the relation (j)n-m{'r, 0) — Note that a characteristic length Iq — \Jh/mVt sets the range 

of these wave functions. 



APPENDIX B: EXCHANGE COUPLING J IN THE HEITLER-LONDON APPROXIMATION 



Here we write down our analytical results in our calculation of the exchange coupling J in the Heitler-London 
approximation. B and Bdif define the inhomogeneous external magnetic field in the system (Fig. 1). As described in 
section HI A, J is calculated by subtracting the average energies of the Heitler-London states '^s/t (Eq. (11)) leading 
to 

J = {l^l' [(*i?|/^l*fl> + {^L\h\^L) +D]- [2Re{S{^nmL)) + /]} , (Bl) 

where h is the single particle Hamiltonian (3), and S — (^'lI^t?) is the overlap. We also define the direct (D) and 
exchange (I) Coulomb integrals 

D = I ""^'^ , |vl/^(2)|^rfxidx2, (B2) 

J Xl — X2 
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:*L(2)*fl(l)dXidx2, 



(B3) 



where e is the electron charge and e the dielectric constant of the material. We list below all the analytic expressions 
for those quantities as a function of the effective mass m, the characteristic lengths in the right and left wells Iqb./l — 

yjTi/m^^'l^ ^ the magnetic fields Bj^ji^ = BzLB^if, and the confining potential parameters a, Ix, ly, hx, hy (Eq. (H))- 
Note that am is the minimum location of each well, which is in general different than the parameter a. Throughout 



the expressions we use the zeroth order modified Bessel function Jo (a;), 7 
the overlap S, which are listed below. 
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For the average of the single particle Hamiltonian h in the right side, just substitute L for R and vice-versa in the 
expression above. 
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A particularly interesting expression is the homogeneous field exchange coupling J{B, B^if — 0) ||T^, which we list 
below (in addition to Eq. (p^)) which is the inhomogeneous case J{B = 0,Bdif))- Since for this case Br = Bl, we 
simplify our notation by using 7o — Ior/l for the characteristic length. Also the overlap becomes S{B,Bijif = 0) = 
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Ba,f [Tesla] 


■e^ \{00\ AH \nm)\'^ 

n— 1,2,3 m 


1(00| Ai? |nm)|^ 

71 — 2,3 m 


Y^Y^ 1(001 A//|nm)|" 
2^2^ IBoo-^nml 

n — 3 m 





meV 


meV 


meV 


1 


0.041 


0.037 


0.002 


2 


0.163 


0.143 


0.011 


4 


0.744 


0.591 


0.053 


6 


1.95 


1.37 


0.156 


B = 4 Tesla 





meV 


moV 


meV 


1 


0.038 


0.034 


0.002 


2 


0.153 


0.141 


0.009 


4 


0.694 


0.603 


0.041 


6 


1.91 


1.49 


0.142 



TABLE I. Energy correction of the lowest Fock-Darwin state (S) when P, D and F levels are included in an inhomogeneous 
field calculation. AH contains all the inhomogeneous field terms that show up in the one particle Hamiltonian hi (Eq. jl^)). 
The quantum dot parameters (defined below Eq. (^)) are the same as the ones used in Figs. 2, 3 and 4. 



n 


m 


Eigonstatc 










1 


1 


^ . 


2 







2 


2 





TABLE II. S, P and D Fock-Darwin energy eigenstates for positive magnetic quantum number m. For negative m, complex 
conjugate the corresponding wave function. 
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FIG. 1. Confining potential (defined in Eq. (g) of the text) that models the double dot structure we study. The magnetic 
field profile is such that the difference between the fields on the right and left is 2Bdif and the average B = {Br + Bl)/2 is 
in the midpoint between the dots. The quantum dot parameters for all the calculations in this paper are given in the text 
following Eq. (||). 
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B [T] 

FIG. 2. Exchange energy J of the coupled dots in a homogeneous magnetic field B. The molecular orbital approximations 
(S, SP, SPD) are very close to each other and represent a quantitative improvement over the simpler Heitler-London (HL). 
Note that J is plotted in the unit of meV, while the overlap S is dimensionless. The analytical formula for the HL theory are 
given in Appendix B. 
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FIG. 3. Exchange energy J of the double dot structure when an inhomogeneous magnetic field is present with average 
B = (Br + Bl)/2 equal to (Fig. 3a) and 4 (Fig. 3b) Tesla respectively. S^j/ = {Br — Bl)/2 is half the difference between 
the fields at the center of each dot. 
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FIG. 4. Exchange energy J of the double dot when an inhomogeneous magnetic field is present with average B equal to 
(Fig. 4a) and 4 (Fig. 4b) Tesla respectively. For completeness, we show J as a function of Bdif up to 25 T, which is much 
higher than what can be obtained in the laboratory. 
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